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We propose an explanation of the anomalous growth of plasticity in ferromagnets 
near the Curie point. We demonstrate that this effect is caused by spin-dependent 
detachment of dislocations from obstacles under an influence of the internal mag- 
netic field. Magnetization fluctuations grow in the vicinity of the Curie point, yield 
an increase of the detachment probability and, hence, an increase of the plasticity. 
We apply this model for a description of the temperature behaviour of the criti- 
cal stress in nickel and of the microhardness of gadolinium. An external magnetic 
field suppresses the magnetization fluctuations and, hence may suppress the above 
singularities. 

PACS numbers: 

I. INTRODUCTION 

Various experimental data indicate a strong influence of the magnetic structure of ferro- 
magnets on their plastic properties. Transition to a magnetically ordered state is accompa- 
nied by a strong change of plasticity (Zackay and Hazlett, 1953, Nabutovskaya, 1969,1971, 
Wolfenden, 1978, Boiling and Richman, 1969, Echigoya etal, 1973, Echigoya and Hayashi, 
1979, Flor etal, 1980, Retat etal, 1985, Retat, 1987, Maksimova and Maiboroda, 1992, Gu- 
layev and Svistunova, 1996). Measurements with small temperature intervals revealed an 
increase of plasticity in a close vicinity of the Curie point (Zackay and Hazlett, 1953, Nabu- 
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tovskaya, 1969,1971, Wolfenden, 1978). This fact by itself is not so astonishing since many 
properties of ferromagnets change near the Curie point (Vonsovskii, 1974), however, the gi- 
ant value of the observed effect is really surprising. According to Nabutovskaya (1969,1971) 
the microhardness of gadolinium changes near Tc by a factor of two, whereas its elastic 
constants vary by a percent or so. For example, the Young modulus of gadolinium decreases 
in the Curie point only by 1.6% (Spichkin etal, 1999). 

Recently we considered a mechanism of the influence of a magnetic field on plasticity of 
some nonmagnetic crystals (Molotskii and Fleurov, 1997, Molotskii, 2000, and references 
therein). The magnetic field induces transitions (singlet to triplet) between different spin 
states of the radical pairs formed by dangling bonds of a dislocation cores and paramagnetic 
obstacles in the course of the pair formation. The state of the pair may be either bonding 
or antibonding, depending on its spin configuration. The magnetic field can influence the 
relative occupations of these states and lead to an increase of probability of the dislocation 
detachment from paramagnetic obstacles. The crystal plasticity, as a result, increases. 

This approach appeared to be very successful and allowed for an explanation of a large 
number of plasticity related phenomena in nonmagnetic crystals. We propose here to ap- 
ply the same approach when discussing plasticity of ferromagnets. In this case, the role of 
the external magnetic field can be played by the local magnetic induction created by the 
spontaneously magnetised surrounding atoms. Recently, we demonstrated that this internal 
magnetic field may be important for the Invar hardening, leading to an increase of the criti- 
cal resolved shear stress of Invar alloys with a lowering temperature (Molotskii and Fleurov, 
2001). We expect to observe a temperature variation of various plasticity related charac- 
teristics of the crystal connected with the variation of the spontaneous magnetization with 
the temperature. When approaching the Curie point the role of the magnetization fluctu- 
ations may become more pronounced. The typical time of relaxation of the magnetization 
fluctuations is usually much larger than the characteristic times of formation or rupture of 
a dislocation - obstacle bond. Therefore we should consider this process under the influence 
of slowly varying magnetic inductions due to fluctuating magnetization. We will show below 
that this may result in a strongly increasing plasticity in the vicinity of the Curie point. 

It is worth mentioning here that the theory, we develop in this paper, will not take into 
account an interaction of dislocations with domain walls, which may serve as efficient pinning 
centres for dislocation (Seeger etal, 1964). This assumption is justified as long as the typical 
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domain size essentially exceeds the dislocation path length. According to the experiments 
of Zackay and Hazlett (1953), and of Nabutovskaya (1969, 1971) this path is hmited by the 
average distance, L/ ~ 1/ ^/pj, between the forest dislocations. The density of the forest 

dislocations in those experiments was p/ ~ lO^cm^^, hence the dislocation path length could 
not exceed L/ ~ 3 x 10~^cm, which was two to three orders of magnitude smaller than the 
typical domain sizes. 

II. INFLUENCE OF INTERNAL MAGNETIC INDUCTION ON PLASTICITY 

OF MAGNETS 

Internal magnetic field of a ferromagnet due to spontaneous magnetization is rather large 
(~ IT) and is capable of influencing plastic properties of crystals. Such an important 
characteristic of crystal plasticity as the dislocation path length increases with the magnetic 
induction B as (Molotskii, 2000) 

L{B) = Lo{l + ^) (1) 

where Bq is a constant characterizing the dislocation - obstacle bond. Usually its value lies in 
the range from 0.2 to IT. In the case of nonmagnetic crystals B is the value of the magnetic 
induction created by an external source. In the case of a ferromagnet we do not need an 
external magnetic field since a strong enough internal magnetic induction is always present. 
This internal magnetic induction decreases when temperature approaches the Curie point 
Tc and starts strongly fluctuating. One may expect that in the critical region, at T ^ Tq, 
these strong fluctuations of the internal magnetic induction will strongly influence plastic 
properties of the ferromagnet. 

The magnetic field infiuences kinetics of the formation of the dislocation - obstacle bonds. 
The characteristic times of these processes are typically on the order of 10~^s (Molotskii and 
Fleurov, 1997), which is shorter than the typical times of the large scale magnetization 
fluctuations determining the local internal magnetic induction Bi^t at each bond (see, e.g.. 
Ma, 1976). This allows one to consider the induction Bint at each moment of time as created 
by a particular space distribution of the magnetization fluctuations in the ferromagnet. 

For example, we estimate here the average dislocation free path length. According to 
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equation (|l|) it is determined by the mean square of the internal magnetic induction, 

{BL) = bI, + {AB'). (2) 

Here Bint is average value of the internal magnetic induction, whereas AB = Bint — Bint- 

Fluctuations of the local internal magnetic induction acting on a particular dislocation - 
obstacle bond can be described using the theory of the second order phase transitions (see, 
e.g.. Landau and Lifshitz, 1980). The local internal magnetic indiction of a homogeneously 
magnetised crystal can be estimated with the help of the Lorentz formula 

An 

B = -M. (3) 

The probability that the magnetization M deviates by AM from its average value M in the 
volume V around the bond is 



1 v^*i/2 I expj- ^^^ l, V>V* 

0, V <V* 



2 {2'KkBTxY''^ 



where x is the magnetic susceptibility, is Boltzmann constant. The probability density 
function (D is truncated at volumes smaller than certain characteristic volume V* which will 
later serve as a fitting parameter. The general approach to the order parameter fluctuations, 
applied here, considers the ferromagnet as continuus medium. Hence, it holds at distances 
which are larger than typical interatomic distances. Therefore one may expect the typical 
scale of the excluded volume V* not to exceed essentially several lattice spacings. Similar 
approach is implied by equation @, which is obtained by considering a small empty volume 
surrounded by a magnetised medium. 

The mean square fluctuations of the magnetization can be then calculated as 

(AM^) = £AM{AMyw{AM, V) = (5) 

Now using equations (0) and the mean square local magnetic field acting on the dislo- 
cation - obstacle bond can be obtained using the following equations 

{B,nt = B,nt ' + {AB') = -^M 2 + g ^ . (6) 

An inhomogeneous spatial distribution of the magnetization, caused by fluctuations, re- 
sults in a demagnetization field, which should have been accounted for in equation (|^). The 
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demagnetization field is the strongest for the spherical fluctuations and is negligible for a 
fluctuation having the form of a plate. The exact account of the demagnetization is a very 
tedious task but finally leads to a numerical factor, smaller than one, in the second term in 
equation (^. This factor is absorbed in the volume V*, which anyhow serves as a fitting 
parameter. 

The temperature dependence of the mean square local magnetic field in the vicinity of 
the Curie point can be readily calculated. For this we may use the temperature dependence 
of the magnetic susceptibility which has the form 



X(T) 
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at T > Tc. 



Here 7 is a critical index. 
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is the Curie constant, in which n is the particle density, p is the effective number of the Bohr 
magnetons fiB per atom. 

The first term in equation (^) contains the average spontaneous magnetization M{T) 
whose temperature dependence is 



M(T) 



Mn 1- 
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at T < Tc, 



(9) 



at T >T( 



c- 



Here Mq is the limiting value of the magnetization at low temperature, f3 is the corresponding 
critical index 



III. TEMPERATURE DEPENDENCE OF THE CRITICAL STRESS: NICKEL 

AS A TEST CASE 

An influence of the magnetization on plasticity was first discovered half a century ago by 
Zackay and Hazlett (1953), who studied the temperature dependence of the critical stress 
for nickel. They have found that the dependence of o"c(T) passes a minimum near the Curie 
point Tc- Samples with 99.95% nickel content were used. There was also a small amount 
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of iron (0.03%) and magnesium (0.02%) atoms. The paramagnetic iron atoms are efficient 
obstacles for dislocations in nickel, and they are most probably responsible for the observed 
sensitivity of the plasticity of nickel to the magnetic field. 

The temperature dependence of the critical stress in a ferromagnet can be found, ac- 
counting for the part played by the internal magnetic fields. Dislocations can be bound to 
paramagnetic obstacles with the binding energy Wm- Then the critical stress is (Friedel, 
1964) 

\Wm\ 



hn{B) 



1 



Tn 



(10) 



where h is the value of the dislocation Burgers vector, and / is the average distance between 
the obstacles. The temperature dependent factor in equation ( p!OD accounts for the thermal 
activation processes (Haasen, 1983). The parameter Tq is proportional to the dislocation - 
obstacle binding energy. 

It was shown in our paper (Molotskii and Fleurov, 1997) that the average length of the 
dislocation free segment depends on the magnetic field similarly to (0), 



1{B) = Zo 



(11) 



Considering a ferromagnet we should introduce the fluctuating internal local magnetic in- 
duction Bj„( = Bint + which acts on each particular dislocation - obstacle bond. It 
means that we have to average equation (|ll]) over various values and orientations of the vec- 
tor AB. This averaging for the case of an easy axis ferromagnet is carried out in Appendix), 
so one gets 

T\ 3 



^^^^^ - 7(5, {AB')) 



0"cO 



Tn 



(12) 



The critical stress appears now to be a function of the average spontaneous magnetization 
and of its fluctuations. The temperature dependence of both these quantities is determined 
by equations (^, (^, and (|^). As a result we now have an equation for the temperature 
dependence of the critical stress in the vicinity of the Curie temperature. 

Figure 1 presents a comparison of the temperature dependence of the critical stress in Ni 
calculated by means of equation (|1^) with the available experimental data. The theoretical 
curve is plotted using the following parameters for Ni: Tc = 627K, Mq = 510G, p = 0.606 
(Kittel, 1986), /3 = 0.33, 7 = 1.33 (Kadanoff etal, 1967). The effective number, p, of Bohr 
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FIG. 1: Temperature dependence of the critical stress in Ni. 

magnetons in Ni is rather small and, hence, one gets a small Curie constant C = 7 x 10~^K. 
The remaining parameters. To = 2070K, Bq = 0.49T, V* = 2.64 x lO^^^cm^, are determined 
by fitting equation (|l^) to the experimental data by Zackay and Hazlett (1953) for their 
samples with 2% deformation. One observes a reasonable general agreement between the 
theory and experiment. However, the measurements were carried out with a ~ 20K interval 
between the points. As a result the sharp minimum in the temperature dependence of the 
critical stress near the Curie point, predicted by our theory, might have been overlooked. A 
relatively small, ~ 15%, decrease of the critical stress in this region is indicative of a possible 
existence of a sharper minimum in this temperature range. 
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IV. TEMPERATURE DEPENDENCE OF THE MICROHARDNESS: 
GADOLINIUM AS A TEST CASE 

This section discusses the behavior of crystal microhardness near the Curie point. The 
results are compared with the available experimental data on gadolinium. Currently there 
is no microscopic theory of microhardness, and hence, no microscopic theory of the influence 
of a magnetic fleld on it. Nevertheless, the dependence of the microhardness on a magnetic 
field can be estimated by means of the following simple considerations. It is known that the 
microhardness H varies inversely proportionally to the plasticity — the higher the plasticity, 
the lower the hardness. On the other hand, the plasticity is directly connected to the 
dislocation path length. The magnetic field dependence of L{B) is given by equation (p. 
Hence, we assume that the average microhardness of a crystal as 

H = Hoj^ (13) 

where Hq is the microhardness in the absence of a magnetic field. 
Assuming an isotropic distribution of the orientations one gets 

H = HoJ(B,{AB')) (14) 

where the function f{B, {AB'^)) is calculated in Appendix. Similarly to the previous sec- 
tion, we use equations (H), (^, and (j^) in order to get equation (JT^) for the temperature 



dependence of the microhardness in the vicinity of the Curie temperature. 

A comparison of this theoretical result with the experimental data available for the micro- 
hardness of gadolinium is presented in Figure 2. The Curie constant C for gadolinium can 
be found using equation (|^) and the values of the relevant parameters: n = 3.02 x lO^^cm"^, 
p = 7.10, Mo = 2010G (Kittel, 1986), which lead to C = 0.33K. The Curie temperature for 
gadolinium is Tc = 292K (Nabutovskaya, 1969). The critical indices are (3 = 0.3265 and 
7 = 1.239 (Aliev etal, 1988). 

The value of the parameter Bq could have been possible to determine from the magnetic 
field dependence of the plasticity. Unfortunately, we do not know about such measurements 
in gadolinium. However, magneto- (Al'shits etal, 1990) and electroplastic (Okazaki etal, 
1979) effects have been measured in zinc and titanium, who, similarly to gadolinium, possess 
hep structures. Fitting these results to the theory (Molotskii, 2000) provides Bq = 0.70T 
for Zn and Bq = 0.94T for Ti. 
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FIG. 2: Temperature dependence of the microhardness of Gd. 

The theoretical curve in Figure 2 uses the values Bq and V* as fitting parameters. The 
obtained parameter Bq = 0.87T for gadolinium appears to be closed to its experimentally 
measured values for Zn and Ti. The characteristic volume is chosen Vq = 1.85 x lO^^'^cm^. 
It corresponds to a sphere with the radius 17A, about four interatomic spacings. The results 
are not very sensitive to the variation of this volume, which can be chosen smaller (up to 
2 X 10~^^cm^) without damaging the overall agreement. 

V. CONCLUSION 

Accounting for the internal magnetic induction of ferromagnets and for its fluctuations 
leads us to the conclusion that various plastic properties of crystals should be singular near 
the Curie point. The critical behavior of the magnetic susceptibility is transferred to the 
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temperature dependence of the critical stress and microhardness. It would be interesting to 
measure other plasticity characteristics, such as dislocation path length L or plastic strain 
rate e near the Curie point in order to look for a possible singular behavior, which follows 
from the theory proposed in this paper. For example, one may expect that carrying out 
measurements similar to those of Zackay and Hazlett, 1953, on the temperature dependence 
of the critical stress in Ni, but with a smaller temperature interval a sharp minimum at 
T = 627K may be observed. 

An external magnetic field is known to smear out the phase transition, to decrease the 
magnetic susceptibility, and to suppress fluctuations near Tc- As shown by Dan'kov etal 
(1998) a 0.5T magnetic field diminishes the magnetic susceptibility near the Curie point in 
gadolinium by an order of magnitude. This should lead to a nearly complete suppression 
of the magnetization fluctuations and, hence, to a suppression of the singular temperature 
dependence of microhardness. As for nickel, a stronger field about 3T (Hischler and Rocker, 
1966) is necessary to suppress the fluctuations, hence correspondingly higher field will be 
necessary to influence the singular behavior of the critical stress in nickel near the Curie 
point. 

We may conclude that an experimental observation of a suppression of the above singu- 
larities in the temperature dependencies of plasticity characteristics in a magnetic field will 
be an important experimental evidence in favour of the theoretical approach to the plasticity 
of ferromagnets presented in this paper. 
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APPENDIX A: AVERAGING FOR AN EASY AXIS CUBIC CRYSTAL 

The averaging is carried out assuming a cubic symmetry. It is certainly relevant to Ni. 
As for the hexagonal Gd, its anisotropy is rather weak and the result obtained here, will be 
also applied for this crystal. We have to average the quantity 

f{B) = (Al) 

-Do 

where B = B + AB, over the distribution (^). Disregarding a possible anisotropy of the 
fluctuations of the magnetic induction in the vicinity of the dislocation - obstacle bonds this 
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averaging can be presented as 



' d cos 9 
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Now integration over Ai? is carried out by parts and the averaging (|A2|) is represented 



as 
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where the notations a = Bo/^{AB^), b = B/^{AB'^), v = V/V% and x = AB/^{AB^) 
are introduced. The identity 
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/ dp exp{—p[a ± i{b + x)]} 
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(A5) 



a±i{b + x) 

is substituted into (|A^). Now integration over x is carried out and a new integration variable 
z = p/ y/v is introduced instead of v. After that, p is substituted for pa and one gets 
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is the probabihty integral. 
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The value of f{B, {AB"^)) can be analytically estimated in the limit of small fluctuations, 



when y (AB"^) <^ Bq. The integral (|A6|) converges for small values of p and the expansion 
~ 2x/v^ can be used. It converges to a rather obvious result 



Bl 



b'- 



(A8) 



For large fluctuations when y {AB"^) ^ Bq the integration ( [A 61 ) includes a broad range 

of large values of p. The probability integral $(a;) rapidly converges to one for x ^ 1. Then 

-B ^/2 

we may exclude the range p < po = ^ ^ whose contribution to the integral is of the 



order of 




(AS2) 



and integrate for p > pq. Then one gets 
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where C = 0.577 is the Euler number. Ei(a;) is the integral exponential function. 
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